One-shot quantum measurement using a hysteretic DC-SQUID 
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We propose a single shot quantum measurement to determine the state of a Josephson charge 
quantum bit (qubit). The qubit is a Cooper pair box (CPB) and the measuring device is a two 
junction superconducting quantum interference device (dc-SQUID). This coupled system exhibits 
a close analogy with a Rydberg atom in a high Q cavity, except that in the present device we 
benefit from the additional feature of escape from the supercurrent state by macroscopic quantum 
tunneling, which provides the final read-out. We test the feasibility of our idea against realistic 
experimental circuit parameters and by analyzing the phase fluctuations of the qubit. 



PACS numbers; 03.67.Lx, 74.50.+r, 85.25.Dq 

A Cooper pair box is a controllable macroscopic two 
level system 0, 0, which is considered as a potential 
qubit in the context of quantum computing ^3] . Coher- 
ent Rabi oscillations have been observed in a CPB, us- 
ing ultrafast pulses Recently, in a superconducting 
device named " quantronium" , coherent oscillations were 
observed using microwave pulses Q . The observed oscil- 
lations lived for almost a microsecond, making supercon- 
ducting circuits promising for realizing quantum gates. 

Despite this progress, quantum measurement on a 
CPB still remains a challenge. The charge read-out cir- 
cuits, like the one in or those using single electron 
tunneling transistors (see, e.g., 0) are far from ideal due 
to decoherence induced by the unavoidable background 
charge fluctuations and by backaction during the mea- 
surement, e.g., in form of shot noise. The former prob- 
lem can be avoided to a large extent if the qubit is op- 
erated at the degeneracy point of the two level system: 
at this optimum point linear coupling to charge fluctu- 
ations is absent. In the quantronium experiment, both 
problems were eliminated successfully by measuring the 
two quantum states of a Cooper pair transistor (CPT) at 
the optimum point by a hysteretic Josephson junction. 
This junction, working in its classical regime, measured 
the persistent current of the two quantum states of the 
CPT. This method cannot, however, be used to measure 
the quantum state of a CPB. Moreover, the measurement 
did not resolve the quantum state in one shot. 

In this letter we propose a new quantum measurement 
procedure based on the entanglement between two quan- 
tum systems: a CPB coupled to a superconducting res- 
onator. The dynamics of this coupled system has been 
theoretically investigated recently 0, S "sl • Here we 
show that at the degeneracy the two quantum states of 
the CPB can be resolved in one shot. The read-out de- 
vice is a current-biased dc-SQUID (the resonator) and 
it is controlled by adiabatic pulses of flux. Very high 
sensitivity and fast measurement can be reached with 
this method. A dc-SQUID is preferred over a Joseph- 



son junction because manipulations using external flux 
do not suffer from time limitations. In a Josephson junc- 
tion, nanosecond pulses of bias current cannot be applied 
because low pass filters are necessary to exclude high- 
frequency noise 10]. 

The circuit of a Cooper pair box coupled to a current- 
biased dc-SQUID was inspired by experiments on a Ry- 
dberg atom in a high Q cavity [Tlj. For a bias current 
/ well below the critical current /c, the CPB and the 
SQUID play the roles of the Rydberg atom and the high 
Q cavity, respectively. For example, Rabi oscillations are 
predicted to occur as a result of spontaneous emission and 
re-absorption by the CPB of a single oscillation quantum 
in the SQUID 6J. The device can be considered as a two 
level system coupled to a harmonic oscillator. However, 
for a bias current very close to, but below the critical 
current, macroscopic quantum tunneling (MQT) in the 
SQUID can occur. This quantum escape phenomenon 
has no equivalent in high Q cavity experiments, and it 
introduces a new element into the dynamics of the sys- 
tem. We will describe how we can benefit from MQT 
in order to do a very fast one-shot quantum measure- 
ment on the CPB. We will discuss the performance of 
this measurement and its back-action on the CPB. 

The superconducting quantum circuit that we consider 
is shown in Fig. ^ It contains three parts which we 
describe in detail below: a Cooper pair box, a hysteretic 
dc-SQUID, and a coupling capacitance between the two. 
The circuit is connected to the external classical circuits 
by three different couplings: a resistance R in parallel to 
a DC current source and voltmeter, a mutual inductance 
to a source of flux modulating current pulses, and a gate 
capacitance to a dc and pulse voltage source. 

The Cooper pair box consists of a small superconduct- 
ing island, coupled to a gate-voltage by a gate ca- 
pacitance Cg. The island is furthermore connected ca- 
pacitively to a superconducting electrode via a Joseph- 
son junction with a capacitance Cj and Josephson en- 
ergy E}. We are interested in the limit of small Joseph- 
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FIG. 1: The CPB-qubit is coupled to a measuring SQUID 
through the capacitance Cc. Vg controlled the CPB states, 
/ and I(t) the measurement procedure. The voltmeter Vm 
performs the read-out. 



son energy, Ej < Ec,j, where Ec,j = 2e'^/Cj^cS is the 
elementary charging energy of the box with Cj,cff = 
Cj + [1/Cs + l/(Cc + Cg)]-i; C5 is the SQUID capaci- 
tance and C'c the coupling capacitance. Let us introduce 
the basis of charge states |n(q)), where n(q) corresponds 
to the number of excess Cooper pairs on the island. If the 
dimensionless gate charge Ng = — CgV^/(2e) is close to 
1/2, the charge states |0(q)) and |l(q)) are almost de- 
generate and the relevant eigenstates of the CPB are 
superpositions of these charge states. Specifically, the 
ground state and the first excited state are, respectively: 

|-> = (|0(q)) + |l(q)))/V2 and 1+) = (|0(q)) - |l(q)))/V2. 

The corresponding eigenenergies are Ezp = i?c,J T ^-E-J- 
We see that the CPB effectively behaves as a quantum- 
mechanical two- level system 's^. 

The hysteretic dc- SQUID consists of a superconducting 
loop with two underdamped Josephson junctions which 
both have ideally the same critical current Iq and ca- 
pacitance Cs/2. We neglect the loop inductance and 
effects due to asymmetry in the SQUID. The damping 
of the two junctions is also neglected now, but it will 
be discussed at the end. With these approximations, 
the SQUID equation of motion is similar to that of a 
single Josephson junction, describing a particle of mass 
TO — Cs,ofr(*i'o/(27r))^ in a tilted washboard potential 
U{ip) = Es[—If/Ic — cos{(f)], where is the phase dif- 
ference of the SQUID, / the bias current through the 
SQUID, $0 = h/2e the superconducting flux quantum, 
<1>DC the external flux, and Ic = 2 |cos(7r<I>Dc/'J'o)| ^0, 
Es = /c*o/(2^) andCs,eff = Cs + [1/Cj + 1/{C, + Cg)]-' 
the effective critical current, Josephson energy and capac- 
itance of the SQUID, respectively. 

For values of bias current not too far below the critical 
current, the potential U{ip) can be well approximated by 
a cubic potential. The quantum dynamics of the SQUID 
using the reduced momentum and position operators 
P = {1/ y/mhujp)P and X — {y^mujp/h)X , respectively, 
is then described by Hq = ^hujp{P^ + X^) + a{I)hujpX^ , 
where X = 1^ is the phase difference, P its conjugate 
operator and cUp = [2^/,/($oCs,eff)]'/'[2(l - I/QV^^ 
the effective plasma frequency of the SQUID. The pa- 
rameter cr(J) = -l/(6a)[2(l - ///c)]"^/"*, where a = 
;l-l/2[$o/(27^)]3/4(Cs,off/c)^/^ gives the relative magni- 
tude of the cubic term as compared to that of the har- 



monic oscillator term. 

For values of / below Ic with a{I) ^ 1, the potential 
barrier is high compared to hup , and the cubic term in Hq 
can be neglected. Many low-lying states are found near 
the minimum of the quadratic potential. The broadening 
of these states due to tunneling can be ignored. Hence 
these states are well approximated by harmonic oscilla- 
tor eigenstates, denoted by |0), |2), corresponding 
to the presence of 0, 1, 2, ... oscillation quanta in the 
SQUID, respectively. Thus, at low enough bias current, 
the SQUID behaves as a superconducting quantum res- 
onator. Since in this limit the phase is localized in a well 
defined minimum of the potential U, the time-averaged 
voltage Vyn across the SQUID remains zero. 

If the bias current is increased such that I ^ Ic, cr{I) 
is no longer negligible and the cubic term affects the 
SQUID dynamics. The barrier height, given by AC/ — 
4\/2/3i?s(l — -f/^c)'^^^, and ujp decrease and vanish at the 
critical current. The number of localized states in a given 
well decreases. Moreover, the remaining energy levels 
broaden due to quantum tunneling from the metastable 
wells of the potential U . The broadening of the ground 
state energy for low damping is given by the tunnel- 
ing rate Fo = a;p6A/6/7r-\/AC//?iajpexp(— 36At//(5fiwp)). 
Since the excited states \n) are located closer to the top 
of the barrier, the tunneling rates and hence the broad- 
ening of these states increase with increasing energy of 
the level as F„ ~ (exp(— 36/5))"ro . Specifically, the 
tunneling rates from states |0) and |1) of the SQUID are 
different by approximately a factor thousand 0, [T3| . 
After a tunneling event, the phase of the SQUID is no 
longer localized and the time-averaged voltage Kn across 
the SQUID is finite. 

The magnetic flux through the SQUID affects its crit- 
ical current Ic, hence Wp, cr(J) and therefore Hq de- 
pend on the flux. For small, time-dependent variations 
of flux, (5<I>(t) ^ $0, and for SQUID parameters such 
that a 3> 1, the total time-dependent Hamiltonian is 
given by H{t) ^ Ho - A(/, ^uc){S^{t) /^o)f^pX , where 

A(/,$Dc) = a[2(l-///c)]-3/8j//,^tan(^$Dc/$o). Be- 
low we consider the effect of a time-dependent perturba- 
tion S^{t) on the dynamics of the coupled system. 

The coupling capacitance Cc plays a crucial role in 
the circuit of Fig. ^ since it couples the qubit and 
the SQUID to each other. Physically, Cc couples the 



charge 



A^g on the CPB to the charge on the 



SQUID. The coupling Hamiltonian can be written as 
He = ~iEcoup\in(^q))P. Here we introduced the char- 
acteristic coupling energy i?coupi = ^/f^p/ Ec,sEc,c/4: 
where Ec,c = 2eVCc.off and Ec,s = 2e^/Cs,ce with 
Cceff = Cc+Cg+{'^/Cj+l/Cs)-'[{Cj+Cs)/iCc+Cg)]/2. 
This coupling energy leads to full entanglement between 
the states of the CPB and the SQUID at the resonance 
condition Ej = fuvp |^. 

Having detailed the superconducting quantum circuit, 
we will now describe the measurement procedure. Sup- 
pose that, at time t = 0, the CPB is in a coherent su- 
perposition a|— ) + as a result of quantum oper- 
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FIG. 2: The quantum measurement procedure is illustrated 
for an initially state |+). 



ations performed at times t < using the gate-voltage 
Vg 0- We propose a way to measure the probability 
|/3p. The measurement procedure is depicted schemati- 
cally in Fig. 121 and consists of three successive step-like 
variations of flux through the SQUID. The steps are not 
sharp, but have a finite rise and fall time St. The flux 
steps must be adiabatic in terms of the dynamics of the 
CPB, 5t h/Ej, and of the SQUID, 5t > 2tt /ujp. But 
they must be instantaneous in terms of the coupling dy- 
namics, 6t ^ h/Ecoupi- The first step at t = puts 
the SQUID into resonance with the CPB during a time 
Tq. The second step at Tq drives the SQUID close to its 
critical current during a time At. The last step sets the 
SQUID far below the critical current. 

In more detail, at t < 0, during the quantum manipu- 
lation of the CPB, the SQUID must be decoupled. This 
condition is achieved if {hujp — Ej) 3> fi'coupi, i-e., off- 
resonance. Thus, in leading order, the eigenstate of the 
entire system is a product of the eigenstates of the qubit 
and the SQUID, in spite of the presence of the coupling 
capacitance. At < = 0, suppose the SQUID in its ground 
state |0). The quantum state of the entire system is then 
|V'(t = 0)) = [a|-)+/3|+)]®|0). 

The flux step applied at t — reduces the effective 
critical current and the resonant condition Ej = fkUp is 
satisfied for a time Tq. The state |— ) fX" |0) is still sta- 
tionary at this resonance since \Ej + ftwpj 3> -Ecoupi- But 
|-|-) ig) |0) is no longer an eigenstate: the state of the 
coupled system oscillates in time between |-|-) ® |0) and 
|— ) (g) |1) at the angular frequency 2£'coupi/S.. Thus after 
a time Tq = /i/(4i?coupi), |+) ® |0) has been transformed 
into |— ) (g) At this point, a second flux step is ap- 
plied through the SQUID which drives the system out 
of resonance. The dynamics is therefore " frozen" in the 
superposition = Tq)) = |— ) (E) (ajO) +e*'''/?|l)) where 
rj is the relative phase arising from the evolution of the 
initial qubit state during time Tq. The full entanglement 
has transferred the coherent superposition of the CPB 
to the SQUID, i.e., the information on the initial qubit 
state is now contained in the SQUID. 

The second flux step reduces the effective critical cur- 
rent such that the constant bias current is close to Ic. 



The barrier is therefore significantly decreased and the 
tunneling rates Fq and Fi are drastically increased. Dur- 
ing the time At, the SQUID can escape from the well to 
the non-zero voltage state by tunneling. If At satisfies 
l/Fi <ti At l/Po, and relaxation between the levels is 
neglected, the SQUID is in its finite-voltage state if and 
only if the SQUID was in the state |1) at time Tq. In 
other words, the proposed measurement determines the 
state of the SQUID in one shot. The escape probability 
corresponds to the amplitude of the initial superpo- 
sition in the CPB. The lack of perfect contrast between 
the escape rates from the two states and relaxation pro- 
cesses introduce an intrinsic error in the proposed quan- 
tum measurement procedure; the influence of this will be 
estimated later. 

At t = Tq + At, the flux is switched back to its ini- 
tial value in order to prevent further tunneling. Because 
the dc-SQUID is hysteretic, the zero-voltage and finite 
voltage states are stable for sufficiently long a time to 
perform the read-out. The first two steps of duration 
To and At perform the quantum measurement. The last 
step provides the classical read-out measurement. 

To check the feasibility of our measurement procedure, 
we use some typical values for parameters of an alu- 
minium superconducting circuit. For the CPB, we choose 
Ej = 26.2 fieV, Cg = 10 aF and Cj = 0.63 fF, and for 
the SQUID, Iq ==°1 fiA, I = 1.1 fiA, Cs = 1 pF and 
$Dc/$o = 0.277. For 5^t))/^Q = 0.013 during the 
first step, the resonance condition is satisfied with 5 lev- 
els in the well; the escape time from the level |1) is much 
longer than 1 ms. An additional increase of flux (5$(t)/$o 
by the same amount is enough to change hiOp by 4 /ieV 
(3> i^coupi) and the system is driven out of resonance. 
The escape time of level |1) drops to about 1 ns. Finally 
we choose Cc = 0.1 fF yielding i?coupi = 0.21 /xeV. Us- 
ing the measuring times Tq = 4.5 ns and At = 5 ns, the 
one-shot quantum measurement can be performed. 

We now turn to the influence of relaxation in the 
SQUID on the measurement procedure. At low temper- 
ature the rate Fr of relaxation between adjacent levels 
down due to interaction with the environment is domi- 
nant. Assuming Fi ^ Fo,F/j, the escape probability at 
t = Tq + At is given in the lowest order by 

^a|0>+;3|l> = l/^l' + l"l'[l-e"'''^1 

-|/3p[e~ri^* + (FR/Fi)(e-roAt _ g-r.At)]^ 

Neglecting the influence of relaxation in Eq. and as- 
suming an infinite contrast between Fq and Fi, the es- 
cape probability gives the amplitude of the initial 
superposition in the CPB. The finite contrast between 
the two states introduces an intrinsic error in the pro- 
posed quantum measurement procedure which is about 
0.8%. Taking into account the relaxation processes, the 
escape probability versus the pulse duration At is plotted 
in Fig. fusing the experimental parameters listed above 
for three different initial states : the two states |— ) and 
\+) and the coherent superposition (|— ) -I- \+))/\/2. The 
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FIG. 3: Escape probability for three different initial states : 
h) , 1+) and (h) + l+))/V2 for Q = 500. The parameters 
for the circuit used in the calculation are given in the text. 

relaxation rate is defined as Tn = lo^jQ where the Q- 
factor was chosen to be Q = 500. The measurement of 
the escape during Ai = 5 ns is a direct measurement of 
the |/3p amphtude. For Q = 500, the total error of the 
one-shot measurement is less than 4%. The presented 
values of 0.8% and 4% error are overestimates because 
the semiclassical approximation of the tunnelhng rate 
slightly underestimates tunnelling from the excited state 
in particular. 

It is important to maintain coherence of the qubit. 
Specifically, we will discuss the backaction, i.e., dephas- 
ing of the qubit, due to the measuring circuit formed 
by the SQUID and its electrical environment. To ana- 
lyze this we need to study the amplitude of the fluctu- 
ations of the phase across the CPB. The measurement 
environment consists of a parallel LRC circuit, where 
ig = \/2^^{l — I/Ic)^^^^ is the Josephson inductance 
of the SQUID and R describes the dissipative part of the 
circuit. It is straightforward to obtain the phase fluctu- 
ations [13, ((A(/))2) = {[(j){t) - 0(O)]2), where is the 
phase difference across the small junction in the CPB at 
time t, using the fluctuation-dissipation theorem. These 



phase fluctuations are determined by the real part of the 
impedance seen by the small junction. If the environment 
would be purely dissipative (resistance R) , we would ob- 
tain fluctuations which diverge in time. On the contrary, 
the inductively shunted circuit, realized by the SQUID, 
protects the qubit: asymptotically {t oo), the expec- 
tation value of phase fluctuations levels off to 

Here Rq = h/Ae^ is the resistance quantum, and wc as- 
sumed low temperatures fc^T <^ fi/(isCs,off)- Accord- 
ing to Eq. (PI the inductance, Ls, provides a protec- 
tion against dephasing: the dissipative part of the en- 
vironment does not affect the value of ((A(/))^)oo- Using 
the typical experimental values given above we obtain 
((A(/))2) ~ 0.02, which is much smaller than tt, thereby 
demonstrating the weakness of the residual phase fluctu- 
ations. The low temperature condition is verifled if T <C 
1 K. As we have shown, the SQUID provides protection 
of the qubit from the decoherence induced by the en- 
vironment. Finally, the proposed measurement can be 
realized at the optimum point, where the background 
charge induced decoherence is largely suppressed. 

In summary, we have shown theoretically that a two 
junction SQUID can perform a single shot quantum mea- 
surement of a Josephson charge qubit. We discussed the 
limits of this detector posed by the flnite contrast in mea- 
suring the quantum state of the SQUID, the flnite quality 
factor of the SQUID, and coupling of the qubit to envi- 
ronment noise. 
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